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Abstract 

Coherent diffraction imaging (CDI) of single molecules at atomic resolution is a major goal for 
the x-ray free electron lasers (XFELs). However, during an imaging pulse, the fast laser-induced 
ionization may strongly affect the recorded diffraction pattern of the irradiated sample. The 
radiation tolerance of the imaged molecule should then be investigated 'a priori' with a dedicated 
simulation tool. The continuum approach is a powerful tool for modeling the evolution of irradiated 
large systems consisting of more than a few hundred thousand atoms. However, this method follows 
the evolution of average single-particle densities, and the experimentally recorded intensities reflect 
the spatial two-particle correlations. The information on these correlations is then inherently not 
accessible within the continuum approach. In this paper we analyze this limitation of continuum 
models and discuss the applicability of continuum models for imaging studies. We propose a 
formula to calculate scattered intensities (including both elastic and inelastic scattering) from the 
estimates obtained with a single-particle continuum model. We derive this formula for systems 
under conditions typical for CDI studies with XFELs. 
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I. INTRODUCTION 

Coherent diffraction imaging of single molecules at atomic resolution [1] is a major goal 
or the present and forthcoming x-ray free electron lasers (LCLS, SACLA, European XFEL) 
4|. The high intensity and the extremely short pulse length (< 100 fs) of the XFEL beam 
are needed to get a sufficiently strong scattering signal from the sample and to reduce the 
effect of radiation damage (in comparison to low fluence imaging experiments, where the 
damage accumulates over many shots). The sample is destroyed during each measurement 
(so-called diffract-and-destroy method [l]). A single diffractive pattern is so noisy that 
averaging over patterns corresponding to the same molecular orientation is needed in order 
to increase the signal-to-noise ratio [5|4lQfl. So for a 3D reconstruction, many patterns of the 
same molecule are needed. 

X-ray diffraction gives information on the electron density of the system, so any changes 
of the electron density due to the radiation damage influence the diffraction patterns. There 
are three main damage components: (i) atoms are 'loosing' bound electrons due to ioniza- 
tion, (ii) ionized atoms may move from their original positions, and (iii) scattering on the 
increasing free-electron density contributes to the signal. In order to investigate the effect 
of these damage components, detailed modeling and understanding of ionization dynamics 
are needed. The continuum approach (Boltzmann ll|, hydrodynamic [12|) is an efficient 
way of modeling damage within large systems. However, the presently available continuum 
methods follow the average single-particle densities, and the recorded diffraction intensities 
reflect two-particle correlations. The information on two-particle correlations is inherently 
not accessible from a single-particle continuum approach. 

In this paper we analyze this limitation of continuum models and discuss the applicability 
of continuum models for imaging studies. In order to study this specific problem we do not 
need to construct a full 'ab initio' model of radiation damage. We analyze this problem 
with a simplified molecular dynamics model reproducing the conditions typical for coherent 
diffraction imaging studies that are known from experiments and from the theoretical studies 
of radiation damage [lJ]. The structure of the paper is as follows. In section II we define the 
relation between total (elastic and inelastic) diffracted signal and its corresponding estimate 
from the continuum model. In section III we describe the physical conditions developing 
within XFEL irradiated samples. In section IV we derive an approximate formula for signal 
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scattered off an XFEL irradiated sample. It is derived from the estimates obtained with 
single-particle densities. This formula is then tested against numerical simulations in section 
V. The results obtained are also discussed therein. The Appendix contains information on 
the simulation details. Finally, in section VI our conclusions are listed. 



II. IMAGING STUDIES AND CONTINUUM MODELS 



The state of an imaged sample during a CDI experiment is dynamic, it changes with time 
as a consequence of the sample irradiation. Stochastically occurring ionizations change the 
electronic states of the atoms at different times, the density of the released free electrons 
increases. Therefore the results of a CDI experiment, even when repeated under the same 
conditions, will vary, due to the stochasticity. We call a unique time evolution that occurs 
during one CDI experiment a single realization. 

In a single realization the total (elastic and inelastic) x -ray intensity scattered off a system 



of electrons in state |^) at time t is proportional to 13|, Il4j : 



J( g ,t) = Jd 3 rd 3 r' (V\n(r,t)n(r',t) \V) e iq - {r ~ r ' ] = {9?\n(q,t)n*(q,t) (1) 

where n(r,t) and n(q,t) are the electron density operator, and its Fourier-transform, re- 
spectively: 

n(r,t) = ]T5(r-f,(t)), 
j'=i 

r N el 

n(q,t) = d 3 rn(r,t)e iq - r = ]Te"^ {t) , (2) 
' j=i 

and the operator Vj(t) is the position operator in the Heisenberg picture. 

Assuming the coherence time of the pulse to be short compared to the timescales of 
the processes occurring within the irradiated sample, the time integrated intensity can be 
approximated as an incoherent sum of the intensities scattered at instantaneous snapshots 
of the system. The number of photons scattered at the vector q during a single realization 
is then proportional to: 

J(q) = Jdth(t)I(q,t). (3) 

The function h(t) describes the average temporal envelope of the pulse, i. e. the modulus 
of the pulse amplitude squared, which is ensemble-averaged over XFEL shots. During an 



imaging experiment a large number of patterns from single shots is collected and summed 
up, in order to obtain an average pattern. Therefore the average of the recorded signal over 
the measured realizations (R) must be formed: 

(%)>* = Jdth{t) Jd 3 rd 3 r' (h(r,t)n(r , ,t)) R e i '>< r - rr >. (4) 

For operators the symbol {...)r implies the averaging of the expectation value over realiza- 
tions, i. e. for an observable A: 

(A) R = J2Pn(t)^n\A\^ n ), (5) 

where p n (t) is the probability that the system is in one of the realized states |^ n ). For a 
classical quantity B averaging over realization simplifies to: 

(B) R = "£p n (t)B n (6) 

where the quantity, B, assumes a value, B n , with a probability, p n (t), at time t. Here we 
would like to emphasize that our analysis is not restricted to systems in thermal equilibrium, 
so that the different realizations are not corresponding merely to the thermal fluctuations of 
an equilibrated system but may also lead through the non-equilibrium stages of the system 
evolution. 

The continuum models, used to simulate the evolution of the irradiated systems, describe 
dynamical properties of electrons and ions, using average single-particle densities (h(r,t)) R . 
The intensity that one can construct directly from these average densities is: 

J c (q) = Jdthit) Jd 3 rd 3 r' (n(r,t)) R (n(r ',*)>* e iq < r ~ r '\ (7) 

The difference between (J(qr))_R and 1 c (q) depends on two-particle correlations during indi- 
vidual realizations. 



III. EVOLUTION OF X-RAY IRRADIATED BIOMOLECULES UNDER THE 
CONDITIONS TYPICAL FOR IMAGING STUDIES 

Detailed simulations of radiation damage jl[ Q, 15, 16] reveal the typical ionization 
dynamics within an x-ray irradiated biological sample during an imaging experiment. Bio- 
logical samples are typically built up from light elements like C, N, O, H. Only a few atoms 



of heavy elements are present, e. g. CI, S, Fe. The primary events that initiate radiation 



damage are the photoionization events 



17|. Photoelectrons released by x-rays from light 



elements have high energy (~ lOkeV). They can leave the irradiated system within a few fs, 
so that the net charge of the system starts to increase. During inner-shell photoionization, 
core holes are created. These holes are then filled by electrons from the outer shells within a 
few tens of fs. For light elements Auger decay is dominant [l^] • Auger electrons have lower 
energies (~ few hundred eVs), and they can be trapped within the charged sample. These 
quasi-free electrons cause further ionizations that end up in secondary electron cascading 
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19|. 



It is energetically preferred for the quasi-free electrons to move towards the center of the 
sample, forming a net-neutral core with ions. This leads to the formation of a positively 
charged outer shell of ions around the neutral core. For samples large enough and sufficiently 
highly ionized, even high-energy photoelectrons can be trapped in the core, enhancing the 
temperature of the quasi-free electrons. 

Conditions of the electron plasma can be characterized by comparing the average potential 
energy of the electrostatic interaction between the electrons, Ecouiomb, to their temperature, 



k B T 
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2l|: 

k B T ~ Ane k B T V 3 J ' U 
where n is the density of the electrons. In case of T > 1, Coulomb interaction dominates the 
electron dynamics and the plasma is strongly coupled. If T < 1, the Coulomb interaction is 
suppressed, and the plasma is weakly coupled. 

Another parameter, the ratio of the electron temperature and the Fermi energy, B, gives 
information on the quantum effects within the system 20|, l2l| : 

e = k JL = k JL f9 ) 

E F J^(3vr 2 n) 2 /3' 

At electron temperatures higher than the Fermi energy (0 > 1), a classical description is 
adequate, while at ^ 1, a quantum mechanical treatment is needed. 

According to the aforementioned simulations 12j, the typical electron density at the 
center of a biomolecule is around 10 23 cm -3 , the electron temperature is ~ 20 eV, when 
photoelectrons are not trapped, and > 100 eV when photoelectrons are trapped. For the 
values of temperature, T = 20 eV and T = 100 eV, the plasma parameters are V « 0.5, 
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rs 2.5, and T a* 0.1, 12.5, respectively, so that the plasma of quasi-free electrons can 
be treated as classical and ideal. 

It should be noted that ions do not form a plasma during an imaging experiment. As 
the imaging pulses are short (~tens of femtoseconds) compared to electron-ion equilibration 
time (~few ps), ions remain 'cold' during the pulse. Unscreened ions from the outer shell 
can move radially, due to the repulsive Coulomb forces. However, if the pulse is not longer 



than several fs, those unscreened ions will also not relocate during the pulse [12|, [22j. As 
a result, during the imaging pulse, the initially neutral system evolves into a two-regime 
system consisting of 'frozen' ions embedded in a thermalized free-electron plasma. 



IV. IMPROVED DESCRIPTION OF DIFFRACTIVE SIGNAL OBTAINED 
FROM SINGLE-PARTICLE CONTINUUM MODEL FOR XFEL IRRADIATED 
SAMPLE 



In this section we investigate in detail how the scattered signal averaged over the real- 
izations reflects the changing electron density within the irradiated sample. We follow the 
approach proposed by Chihara in Ref . 23[ , and further applied in many works investigating 
x-ray scattering on plasmas, e. g. However, in our case we will apply it to a system 

that can be far from thermal equilibrium at its initial evolution stages. This is in contrast 
to the standard plasma applications that usually consider systems in thermal equilibrium. 

Following [23(, we separate the total electron density operator, n(r, t), into the the bound 
electron component, n b (r, t), and the unbound electron component. In general, one could 
develop a model for describing the dynamics of all electrons, also those ones that escape 
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12| follow the 



from the sample. However, the currently available continuum models 
dynamics of particles only within a restricted volume. Therefore we separate the unbound 
electron density operator into two components. One component refers to the electrons that 
escaped from the sample, n e (r,t), and the other one to the electrons trapped inside the 
sample, n t (r,t). From Eqs. ([1]), ([5]) we then obtain: 

2\ 



R 



(h b (r, t) + n t (r, t) + n e {r, *)) e iq r d 3 r 

n b (q,t)\ 2 ) R + (\Mq,t)\ 2 ) R + (\h e (q,t)f 

+2 Re [ (n b (q, t) h*(q, t)) R + (n b (q, t) n*{q, t)) R + (n t (q, t) h* e (q, t)) R \ 
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= (I bb (q, t)) R + (I tt {q, t)) R + (I ee {q, t)) R 

+2Re [{I bt {q,t)) R + (I be (q,t)) R + (I te {q,t)) R ) , (10) 

where Iij(q,t) = (ty\ni(q,i) n*(q,t) \^f) for i,j = b,t,e. First we analyze the effect of the 
diagonal terms, (I bb (q,t)) R , (I t t(q,i)) R , (I ee (q,t)) R , and later the effect of the cross terms, 
(ht{q,t)) R , (I be {q,t)) R , (I te {q,t)) R . 



A. Bound electrons 



We can express the bound electron density operator as a convolution of atomic electron 
densities, h b j(r,t) and of the atomic position operators, 5(r — Rj(t)) 28]: 

N a 



3=1 



Its Fourier transform then reads: 



h b {r,t) = J2 J d 3 r'n btj (r',t)5(r-r' - Rj(t)). 



(11) 



Na 



n b {q,t) = / h b (r,t)e^ r d 3 r = £ n b>j {q, t) e^ R ^\ 

J 3=1 



(12) 



where A^ a is the total number of atoms, Rj{t) is the position operator of atom j, and 
h bi j(q,t) is the Fourier transform of the density h b) j(r,t). The scattered intensity averaged 
over realizations is: 

(I bb (q,t)) R = (\n b (q,t)\ 2 ) R = £ (^.( 9 ,t)^(g,t)e^^W-^W)) fl . (13) 

j,k=l 

Radiation damage can affect the intensity pattern in two ways: (i) form factors of atoms 
can decrease due to the progressing ionization, (ii) positions of ions can change due to 
the Coulomb interaction within the ionized system. The latter effect can be overcome by 
applying pulses short enough (< 10 fs 22]), which are available at XFELs |29(] . However, 
the reduction of the form factors cannot be eliminated. 

Assuming static atomic positions, (Rj(t)) R = Rj in all realizations, and no cor- 
relation between ionization events of different atoms, so that (n bt j(q,t)fil k (q,t)) R = 
(n b j(q,t)) R (n bk (q,t)) R for j ^ k, the corresponding averaged intensity is: 



Na N a 

(I^(q,t)) R = J2(\n bij (q,t)\ 2 ) R + £ (n bij (q, t)) R (n£ fc (q, t)) R e 

3=1 37^=1 



iq-(Rj-Rk) 



(14) 
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N a N a 

E + E ( <l/ite *)l 2 >« - 1 </i(«, 



j=l J'=l 



3=1 



a (N (t)\ a 

E - y. — iSLcj, + E ( (\fj(q,t)\ 2 ) R - 1 (/;(<?, I 2 ) + IM«,0> 



j=l Z 3 j=l 



i2 



where fj(q,t) = (\& | h b j(q,t) | \&) is the atomic form factor, n b (q,t) = Y^f=ifj(q,t) e iqR j 
is the classical bound electron density. Deriving the second row in Eq. (|14|) . we used the 
following relation: 

(MQ)\ 2 ) R = T,Pn(t)(^n\n b AQ)\^n)(^n\ntM)\^n) 

n 

+ (* n \n bJ (q) (1 - |*n)<*«|)n^(g)|*n)) 
= {\fj(q,t)\ 2 ) R +(I^)R. (15) 

If the binding energy of a bound electron is much less than the energy of the incoming 
photon, the intensity Ij nel scattered inelastically on the atom/ion, j, can be approximated as 
being proportional to the number of its bound electrons, N bt j(t), multiplied by the intensity, 
1 r atomic zi scattered inelastically from the neutral element of atomic number, Zf 



I jinel\ R rinel 

\ 7 ]R ^7 aton 



ry aiomic,Zij 

(16) 

According to Eq. (fT4|) (last row), the total averaged scattered intensity from bound elec- 
trons is then the sum of: (i) the inelastic incoherent scattering component, J2j=i (Ij nel )R, 
(ii) elastic incoherent scattering component containing the variance of the atomic form fac- 
tors of the individual ions, and (iii) the density, \(n b (q,t)) R \ 2 that describes coherent elastic 
scattering. Similar terms appear in the final formula of Ref. 23|], except for that term which 



is a consequence of the externally driven ionization dynamics. 



B. Trapped electrons 

We investigate now the contribution of trapped electrons to the imaging pattern. Under 
the conditions typical for CDI studies electrons can be treated as classical particles. In what 
follows we will then replace the electron-density-operator expectation values, h(r, t)\ty), 
with their corresponding classical equivalents, n(r,t). The classical electron density and its 
Fourier transform, defined for N t trapped electrons, located at positions Vj(t) at time t, then 



S 



arc: 



N t + -^\{n t {q,t)) R \ 



N t 

n t (r,t) = J2S(r- rj (t)), 

3=1 

r N t 

n t (q,t) = d 3 rn t (r,t)e i(ir = Y<e iqr i {t) , (17) 

J 3=1 

and the realization-averaged signal scattered from trapped electrons is proportional to: 

Aft 

(I tt (q,t)) R = N t + £ ( e ^M')-^))) R . (18) 

Using the two-particle distribution function {30], one can rewrite it further as: 

(I tt (q,t)) R = AT i + ^^|d\dV(^(r,t)) /? (n t (r',t)) i? ( ? (r,r',t)e i9 ("') 

AW 
N 

+ \d\d\> (n t (r, t)) R (n t (r>, t)) R (g(r, r', t) - l)e*<r^\ (19) 

where g(r,r',t) is the two-particle distribution function at time t. The time-dependence of 
g(r,r',t) is imposed by the non-equilibrium evolution of the considered system. 

In the second row of Eq. (|T9l) the first term is a g-independent shift given by the number 
of trapped electrons, N t . The second term depends on the average electron density. The 
third term contains the statistical information about the two-particle correlations. This 
information is not accesible within single-particle continuum approach, as considered here. 
The first term and the third one correspond to the incoherent scattering, whereas the second 
one corresponds to coherent scattering on free electrons as in Chihara's formula in Ref. (23 ]. 

Neglecting all correlations in f lT9|) . we arrive at the following estimate: 

(Iu nC ° rr (q,t)) R =N t + ^^\(n t (q,t)} R \ 2 . (20) 

In the above derivations we assumed that the number of trapped electrons, Nt, is identical 
for all realizations at any time t. As the N t distribution results from a stochastic process, 
this assumption is generally not valid. In the general case, we should replace N t with its 
average over realizations at time t, (N t (t)) R . This is an accurate approximation if the 
number of particles in individual realizations is large. For a small number of particles this 
may introduce a large relative error. However, in our case, the contribution of trapped 
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electrons to the overall scattered intensity is small, so in any case it should not introduce a 
significant error. 

We note that the approximation, Eq. (|2"Uj) . could be further improved by taking into ac- 
count correlations. In general, especially in case of short pulses, the electrons may not reach 
thermal equilibrium during the pulse. This would require a time-dependent parametrization 
of g(r, r', t). If the pulse is long enough compared to the electron thermalization timescale, 
the two-particle distribution function for electrons could be approximated by the equilibrium 
one. 

C. Escaped electrons 

Energetic photoelectrons leave the sample soon after the exposure to the x-ray pulse 
starts. However, they may stay within the focus of the beam during the imaging pulse and 
also contribute to the diffractive patterns. 

In order to describe their contribution to the intensity, the same approximate formula as 
for the trapped electrons, Eq. (I20l . can be applied: 



In case of escaped electrons, their density is dilute, e.g., the distances between escaped 
electrons are on average large. Their contribution will then show up mainly at small values 
of q, limited by the inverse of the size of the escaped electron cloud. We explain here that 
the relevant q range for imaging studies is defined by 1/L < q/2ir < 2/d, where L is the size 
of the object (here the size of the electron cloud), and d is the desired resolution. Therefore, 
if the continuum model does not provide information on (n e (q,t)) we can approximate 
the effect of the escaped electrons as: 



which is accurate enough in the high q regime, limited by the inverse resolution, 2/d. 
D. Cross terms 

The escaping high-energy electrons are typically located at large distances from other 
particles, so the potential energy of their individual interaction with other particles is much 




(21) 




(22) 
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smaller than their kinetic energy. Therefore, the Coulomb correlations between bound elec- 
trons and escaped electrons, and between trapped electrons and escaped electrons, can be 
neglected, and the respective cross terms factorize to: 

(I be {q,t)) R = (n b (q,t)n*(q,t)) R = (n b (q,t)) R (n* e (q,t)) R , 

(I te (q, t)) R = (n t (q, t) n*( 9 , t)) R = (n t (q, t)) R (n* e {q, t)) R . (23) 

Following our discussion from the previous subsection, we can neglect these terms in the 
considered high q regime. 

The correlation between bound and trapped electrons, originating from the same initial 
process of ionization, strongly decreases with time as the number of particles increases and 
the particle cloud evolves. However, there are still Coulomb correlations between these 
particles. If we compare the potential energy of the electrons trapped in the field of the ions 



to their kinetic energy at the typical simulation parameters 12J, we obtain: 



T ei sa 0.93 at T = 20 eV, 

T ei « 0.18 at T = 100 eV. (24) 

If the electron-ion system is weakly coupled (T e i < 1), we can then neglect the correlations 
between the different density terms, and the cross term can be factorized as: 

(I bt (q, t)) R = (n b (q, t) n* t {q, t)) = (n b (q, t)) R t)) R . (25) 

Again, the cross terms then depend solely on average densities, which are known from the 



continuum model. A similar cross term appears in Chihara's formula in Ref. 23J] 



E. Intensities from continuum model 



Here we collect the terms that we derived in the previous subsections. The estimate for 
the total scattered signal in the uncorrelated case in the region of high q values is: 

(i uncorr ( q ,t)) R = E ^# g iS±c,z J + Ef(l/i(g.0l a 



+ N t + 



i=i 



1 , 



(fj(Q,t)) R )+\(n b (q,t)) R \ 



N t 

N a 



(n t (q,t)) R \ z + N e +2Re [(n b (q,t)) a (r$(q,t)) R \ 



I C (g,t) + E( {\f^t)\ 2 ) R - (ffat)) 



' R 
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+ 



V ( N b,j(t)}R jine i J_ w7/ ,s , at , iu 

2^ 7 1 atomic,Z 1 TF 1 t Wi 1 ) + iV < + iV e, 



3=1 



N t 



(26) 



where J c is the intensity obtained from the total average electron density within the con- 
tinuum model (cf. Eq. (J7])): 



I c (q,t) = \(n b (q,t)) R + (n t (q,t)) R \ z 

= I°(q,t) + I?{q,t) + 2 Re [(n b (q,t)) R (n* t (q,t)) 



Rl » 



(27) 



and Ifr t (q,t) = (nb,t{q,t)) R . If the correlations within the electronic system are known, 
Eq. (1261) can be improved by including Eq. (I19p . Otherwise, the uncorrelated approximation, 
Eq. ([20]) , should be used. 



V. COMPARISON TO NUMERICAL SIMULATIONS 

In order to investigate the applicability regime of Eq. (126]) . we performed dedicated molec- 
ular dynamics (MD) simulations of electron-ion systems at the conditions corresponding to 
those that develop in the center of irradiated biological samples during an imaging experi- 
ment. For technical details please see the Appendix. 

Using MD simulations one can follow individual realizations of the dynamics within an 
irradiated system. This enables a direct comparison of the time-integrated realization- 
averaged diffraction signal to its estimate, Eq. (126]) . derived from single-particle densities. 
It corresponds to the estimate that could be obtained from a continuum model. From the 
results of MD simulations we calculated: (i) the average intensity from different realizations, 
(I{q,t)) R (Eqs. P),®), (ii) the uncorrelated intensity, (I uncorr (q, t)) R (Eq. (BE)), and (iii) 
the correlated intensity, (I corr (q,t))n (Eq. (Tl9l) . for electrons only). We integrated them 
over time to obtain (3(q)) R , (3 uncorr (q)) r, and {0 corr (q)) R respectively. We also calculated 
J c (g), according to Eq. (jTJ). We investigated the following cases. 

• Scattering from bound electrons only 

We modelled a cluster consisting of 100 initially neutral carbon atoms. The atoms were 
located randomly within a sphere of radius Rq = 7.1 A. Their density was 1/15 A- 3 , 
which corresponds to a typical protein density. The nearest neighbor distances were 
always larger than 1.5 A. The atoms were not moving. The time dependence of the 
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average number of bound electrons per ion was approximated by an exponential func- 
tion in such a way that the average number of bound electrons per atom decreased 
from 6 to 4 during the simulation time. 

Scattering from trapped electrons only 

Electrons were treated as classical particles with regularized Coulomb interaction be- 
tween them. The granular ions were not included. However, an average positive field 
(radial harmonic potential) was added to keep the electrons together. The strength of 
the harmonic potential was determined by the average ion density and the average ion 
charge. The electron dynamics was followed by integrating the equations of motion 
of the electrons. The parameters we used corresponded to the results of the damage 



modeling for biosamples 12j : n atom = 1/15 A -3 (average atomic density), Q ion = 2 
(average ion charge). The number of electrons was N = 200. The electrons filled a 
sphere of R = 7.1 A radius at the density of 2 ■ n atom . We followed 100 realizations 
of the system during T pu i se = 10 fs, from which we calculated the realization average. 
We performed two sets of calculations, one with the electron temperature of 20 eV 
another one with T e = 100 eV, which corresponded to the cases when photoelectrons 
could escape or were trapped within the system, respectively. 

• Scattering from bound and trapped electrons 

This case was just the merging of the above two cases, with the same values of the 
simulation parameters. A regularized Coulomb potential described the interaction 
between granular ions and electrons. The external harmonic potential was not applied 
here. 

Below we discuss the simulation results that we have obtained. 



A. Bound electrons 

The scattered intensity depends on stochastic properties of the ionization process. Fig. 
Q] shows the fraction of ions with different charges created within the sample. Although 
we used a simplified model for describing the ionization (see Appendix), the trends in the 



charge distributions are in agreement with Ref. 



121. 
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FIG. 1: (Color online) Average ion-charge-state populations within the imaged system as a function 
of time. 

In Fig. [2] we show the intensity distribution for scattering from bound electrons along 
a randomly chosen axis in reciprocal space. The (j™ co ™ "(q)) R gives a good estimate of 
(Jbb(g))ij. There is a disagreement between 0^(q) and (Jm>(<?)).r, as there is no inelastic 
scattering included in 0^'{q). The shift between these two intensities is proportional to the 
number of bound electrons at high q. 

Therefore, one should rather use here (^ corr (q r ))_R as an estimate of the scattered inten- 
sity from bound electrons. 

B. Trapped electrons 

The electronic system considered here is in thermal equilibrium. The inverse Debye- 
length: 

k ° - /IS (28 » 

is ~ 1 A ~ l at the electron temperature 20 eV. In Fig. [31 we show the time integrated scattered 
intensity, (3u{q))r, along an axis in reciprocal space, averaged over realizations, as compared 
to the intensity, Jf (q), obtained from the average density, and to the intensity, (J™ corr (g))_R, 
estimated in the uncorrelated case. The dominant correction (to Jf'(q)) originates from the 
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FIG. 2: (Color online) Time integrated intensity from diffractive scattering only by bound electrons 
as a function of the momentum transfer, q, along a random axis in reciprocal space. We show (i) the 
intensity constructed from average single particle densities 3^ , (ii) the scattered intensity {3bb)n, 
and (iii) the estimated uncorrelated intensity (J^ lcorr )_R. 

granularity of electrons, and is given by the total number of the trapped electrons in the 
system (cf. Eq. ( l20i) ). At larger momentum transfers, q > kn, this correction is accurate 
enough and (J" t ncorr (q))# agrees well with {3tt{q))R- However, at lower momentum transfer, 
it breaks down. The reason for this is that in Eq. (|20|) all correlations are neglected a priori, 
and in our simulated case two-particle correlations are present. They manifest themselves 
at low q. As expected, their effect decreases at higher electron temperatures, T = 100 eV 
(Fig. II . 

For an infinite, homogeneous system, the scattered intensity is connected to the radial 



distribution function by the following expression [30] : 

3(q) =N t (l + n t j d 3 r (g(r) - 1) e^) , (29) 
where n t is the homogeneous density of the electrons. In the case of a Debye Huckel plasma, 



the radial distribution function, g(r), is 



30j: 



Apart from the region of q ^ 1/Rq , where the finite size effects dominate, at T = 20 eV 
C^t (<?))# remains in a good agreement with the intensity profile, 3® H (q), obtained from Eqs. 
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FIG. 3: (Color online) Time integrated intensity from diffractive scattering by trapped electrons 
as a function of the momentum transfer, q. The electron temperature was set to 20 eV, and an 
external harmonic potential was applied. We show: (i) the intensity constructed from average single 
particle densities Jp, (ii) the scattered intensity (Ju)r, (hi) the estimated uncorrelated intensity 
(J" t ncorr )_R, (iv) the estimated correlated intensity (Jj t orr ),R, and (v) the intensity obtained from the 
theory of homogeneous, infinite-size, weakly coupled plasmas H . 

( l2"B]|2"§j) . evaluated with the initial simulation parameters, Rq = 7.1 A, n t = A^/(47ri?o/3). A 
better estimate, (J£° rr (q)) r, can be obtained by evaluating Eq. ( TH)j) with the pair correlation 
function, Eq. (!30|) . It then applies well in the entire g-range. 

At the higher temperature (T = 100 eV), the spread of the electron density is wider 
(Fig. [5]). We have then to evaluate Eq. ( !T9|) and Eqs. ( 128H30]) with parameters that correspond 
to a homogeneous sphere with the average radius estimated to be R ~ 8.2 A. Our findings are 
similar: the intensity fit, (^ t orr (q)) r, obtained from Eq. (TlT)|) is in good agreement with the 
total signal. To sum up, (^ corr (q))R overestimates the full scattered intensity at q < k D , 
however, it still gives an accurate estimate of the trapped-electron background in the region 
of high q, important for high resolution imaging. The estimate, {^tt^iQ^R applies well in 
the entire q region. 
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FIG. 4: (Color online) Time integrated intensity from diffractive scattering by trapped electrons 
as a function of the momentum transfer, q. The electron temperature was set to 100 eV, and an 
external harmonic potential was applied. We show: (i) the intensity constructed from average single 
particle densities Jp , (ii) the scattered intensity (Ju)r, (hi) the estimated uncorrelated intensity 
(J" t ncorr )_R, (iv) the estimated correlated intensity (Jj t orr ),R, and (v) the intensity obtained from the 
theory of homogeneous, infinite-size, weakly coupled plasmas H . 




FIG. 5: (Color online) Average density of an electron system containing 200 electrons (solid lines) 
at temperatures, 20 eV and 100 eV, with an external harmonic potential applied. Dashed lines 
show the radius and density values, used for evaluating Eq. (fT9|) and Eq. (j3DJ). 



17 



C. Bound and trapped electrons 



At last we discuss the case with contributions from both bound and trapped electrons. In 
this simulation electrons were moving in the field of granular ions, therefore we can expect 
some change in the behavior of the average intensity scattered by the trapped electrons 
(Fig. [6]) as compared to that scattered by the trapped electrons in a smooth positive field 
(Fig. [3]). Discrepancies show up as small peaks at the q values corresponding to the locations 
of ion peaks (Fig. [2]). These peaks are a consequence of electron-ion correlations resulting in 
the temporary localization of some electrons near granular ions. Generally, the q dependence 
of the curves in Fig. O is similar to that obtained in the case when a harmonic potential was 
applied (Fig. [3]). 

The intensities scattered by both bound and trapped electrons are shown in Fig. [7J The 
offset between the intensities (J(q))^ and c (q) at high q is due to the trapped electron 
contribution, N t , and the inelastic scattering on the bound electrons. For the considered 
case even the intensity {'J uncorr {q)) r that neglects correlations between the trapped electrons, 
approximates well the full scattering signal. In the calculations presented, N e = 0. 

Let us recall here that the structural information is carried by the elastically scattered 
photons. However, in currently planned experiments on non-periodic single objects, inelastic 
scattering will also contribute to the measured patterns. Coherent XFEL sources have been 
reported to produce radiation with a bandwidth of 0.2 — 1 % [2]. At the photon energy of 
E p h = 12 keV, it corresponds to a bandwidth of AE p h = 20 — 120 eV that is comparable with 
the Compton broadening. As a consequence, even if detecting exactly at the photon energy 
E p h, one collects the elastic scattering signal of incoming photons and inelastic scattering 
signal of incoming photons with the initial energy within the bandwidth. Both the elastic 
and the inelastic scattering then contribute to the measured signal. 

We showed that inelastic scattering on bound electrons can have a significant impact 
on the measured intensities: it contributes to the background that reduces contrast of the 
recorded image. This effect is even more pronounced at larger momentum transfers (^y^eai 
and IbbMeal in Fig. Eli). 



At experiments that have already been done at LC 



was negligible. When the sample is a nanocrystal 



iS, the effect of inelastic scattering 
311 ]. the strong coherent Bragg peaks 



dominate over inelastic scattering. Inelastic scattering is also negligible at low resolution 



18 



300 

in 
250 

_d 200 
i_ 
ro 

>150 

■4-J 

£ ioo 

c 

- 50 







1 








■ 




- s; 






* 

9 










/ (zuncorr \ 
/J? 






i . 


/ acorr \ 








J?" 











0.0 



0.5 



1.0 

q/2ir [A" 1 ] 



1.5 



2.0 



FIG. 6: (Color online) Time integrated intensity from diffractive scattering by trapped electrons as 
a function of the momentum transfer, q. The electrons were moving in the field of positive granular 
ions. The electron temperature was set to 20 eV. We show (i) the intensity constructed from average 
single particle densities 3f , (ii) the scattered intensity (3tt)R, (hi) the estimated uncorrelated 
intensity {Jft 100 "') r, (iv) the estimated correlated intensity (Jj° rr }i?, and (v) the intensity obtained 
from the Debye-Hiickel theory . 

experiments on single objects 



32], when intensity at small scattering vectors is collected. 
However, the above predictions show that the effect of inelastic scattering should be taken 
into account, when planning atomic resolution imaging of non-periodic samples. 

Finally, we show the effect of damage on the recorded total signal (Fig. |8b). As expected, 
progressing damage does not change the positions of intensity peaks that correspond to the 
positions of imaged ions. It only changes the (vertical) positions of intensity minima and 
maxima, which reduces the image contrast. 



VI. SUMMARY 



The continuum approach is an efficient tool for the study of the dynamics of large atomic 
or molecular samples exposed to an intense XFEL pulse. However, the presently developed 
continuum models deliver only information on single-particle densities, and imaging studies 
require information on two-particle correlations. 



19 





1600 




1400 








1200 








1000 


i_ 




ra 


800 


> 




'ld 


600 


c 




ai 
-t-j 


400 


c 






200 



0.0 




2.0 



q/2ir[k ] 



FIG. 7: (Color online) Time integrated intensity from diffractive scattering by bound and trapped 
electrons, as a function of the momentum transfer, q along a random axis in reciprocal space. We 
show (i) the intensity constructed from average single particle densities 3 C , (ii) the full scattered 
intensity (3)r, (hi) the estimated uncorrelated intensity (y uncorr } R ^ anc [ (i v ) the estimated correlated 
intensity, (T t r} R . 



Here we have studied in detail the effect of two-particle correlations on CDI scattering pat- 
terns, obtained from electron-ion systems under conditions similar to those expected during 
XFEL imaging experiments at atomic resolution. We derived an estimate for the scattered 
intensities from the single-particle-density continuum model: Eq. (126 p . and demonstrated 
through numerical simulations that it can describe the scattered signal with a good accu- 
racy. Correlation effects manifest themselves only in the region of low q, together with the 
effects of the finite size of the sample. 

Our results have implications for imaging-oriented studies of radiation damage performed 
with continuum models, as they define the limits of applicability of these models for CDI 
simulations. 



Appendix: Simulation details 

The simulations performed in this work were based on the molecular dynamics approach. 
The particles (ions and electrons) were treated as classical, charged point-like particles. 



20 





1600 




1400 








1200 






X! 


1000 






ra 


800 


> 




'1/1 


600 


c 




ai 
-t-j 


400 


c 






200 




(a) 



0.0 0.5 1.0 

q/2n [A" 1 ] 



1.5 



2.0 



1600 




(b) 



g/27r[A ] 



FIG. 8: (Color online) The effect of inelastic scattering and damage on the scattered intensity. We 
show: (a) the signal from the 'ideal' (undamaged) sample without inelastic scattering Jf^eaZ' and 
with inelastic scattering ideal , (b) the signal from the 'ideal' (undamaged) sample (3)r, and the 
signal from the damaged sample (J)_r, both including the inelastic scattering component. 

• Ionization of atoms 

In order to describe the time dependence of the average number of bound electrons 
per atom, we used an approximate exponential form: 



N(t) = (N - N inf ) eM-wt) + N inf , 



(A.l) 
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where N(t) denotes the number of bound electrons per atom at time t. Initially at 
t = fs, N = 6. The fitting parameters, w and N in f, were adjusted to N in f = 3.98 



and w = 0.456 (fs) 1 , following the results from Ref. 12j. Atomic form factors were 



calculated with the XATOM package 
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34). 



Dynamics of trapped electrons 

The motion of the trapped electrons was followed with the Newton equations solved by 
the numerical Velocity- Verlet algorithm. The interaction between the charged particles 



was described by the regularized Coulomb potential, V(r) = 1/yr 2 + r 2 ,, where the 
cut-off parameter, ro = 0.1 A, was used. The simulation timestep was las . 

Simulating the smooth positive background field for trapped electrons, we applied 
an external radial harmonic potential of the form, V(r) = —Dr 2 /2. Initial electron 
positions corresponded to those ones within an electron system in thermal equilibrium. 

Simulating the case with granular ions, we placed the ions at random positions within 
a sphere of 30 A radius. The ion density was set to 1/15 A -3 , which corresponds to 
a typical atomic density within a protein. The minimal atom-atom distance was kept 
larger than 1.5 A. Ionic charge was set to +2. Let us note that the static electronic 
charge was fully neutralized by ions inside a sphere of 7. 1 A radius. 

Calculation of scattered intensities 

We recorded the total number of the bound electrons at each atom/ion at a time 
step. We did not distinguish between bound electrons at different atomic orbitals. 
Therefore, we approximated the scattering factor of an ion with an average scattering 
factor, averaged over different electronic configurations. The scattering factor of a 
trapped electron was 1. 

In our simulation, we considered diffractive scattering only from ions inside the neu- 
tralizing 7.lA-radius sphere (corresponding to the imaged net-neutral electron-ion 
sample), and not from the ions in the positively charged outer shell. 
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